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We point out that electromagnetic one-way edge modes analogous to quantum Hall edge states, 
originally predicted by Raghu and Haldane in 2D gyroelectric photonic crystals possessing Dirac 
point-derived bandgaps, can appear in more general settings. In particular, we show that the TM 
modes in a gyromagnetic photonic crystal can be formally mapped to electronic wavefunctions in a 
periodic electromagnetic field, so that the only requirement for the existence of one-way edge modes 
is that the Chern number for all bands below a gap is non-zero. In a square-lattice gyromagnetic 
Yttrium- Iron-Garnet photonic crystal operating at microwave frequencies, which lacks Dirac points, 
time-reversal breaking is strong enough that the effect should be easily observable. For realistic 
material parameters, the edge modes occupy a 10% band gap. Numerical simulations of a one-way 
waveguide incorporating this crystal show 100% transmission across strong defects, such as perfect 
conductors several lattice constants wide, larger than the width of the waveguide. 



Photonic crystals — structures with periodicity compa- 
rable to the wavelength of light — possess many unusual 
and technologically important optical properties [H 0, 0| . 
The theory of these devices relies on an analogy between 
Maxwell's equations in a periodic medium and quan- 
tum mechanics with a periodic Hamiltonian, from which 
photonic bandstructures arise in a manner analogous to 
electronic bandstructures in a solid. The existence of 
photonic bands has inspired a variety of applications in 
integrated optics; for instance, photonic crystals made 
of magneto-optical (MO) materials, which break time- 
reversal symmetry (T), can be used to construct non- 
reciprocal optical circuits 0, Raghu and Haldane 
Q have recently extended this analogy in a remarkable 
new direction, by predicting the existence of one-way 
electromagnetic modes similar to the chiral edge states 
found in the integer quantum Hall (QH) effect [7j. These 
modes are confined at the edge of certain 2d MO pho- 
tonic crystals and possess group velocities pointing in 
only one direction, determined by the direction of an 
applied DC magnetic field (the time-reversed versions 
are not eigenstates of the system). Due to the absence 
of back-propagating modes, backscattering is completely 
suppressed. This is potentially important for slow-light 
structures, which are susceptible to backscattering [8|, 
among other applications. However, the requirement 
of locating a "Dirac point" in fc-space, as suggested by 
Raghu and Haldane, led these authors to focus on TE 
modes in triangular lattices and thus to gyroelectric ma- 
terials. In realistic gyroelectric materials, the strength 
of T-breaking, characterized by the ratio of (imaginary) 
off-diagonal to (real) on-diagonal elements of the permit- 
tivity tensor e (the Voigt parameter) is at most ~ 10~ 3 . 
As a result, the bandgap is not robust against disorder, 
and the edge modes scatter easily into bulk modes of 
the crystal, leading to significant radiative loss. To our 
knowledge, these one-way edge modes have never been 
observed experimentally, and to do so one would desire 
a bandgap that is orders of magnitude broader. 
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FIG. 1: (Color) Magneto-optical one-way waveguide, (a) Pro- 
jected band diagram. The dispersion curve (red) for the edge 
modes spans the band gaps for the cladding crystals, (b) 
Steady-state field pattern for E z , the out-of-plane electric field 
component. Blue and red represent positive and negative field 
values. A source, indicated by double arrows, is located at 
the interface between a magnetized (+z) YIG crystal (lower 
half- plane) and an alumina crystal (upper half-plane), and 
operates at a mid-gap frequency (0.552 ■ 2nc/a). The excited 
edge mode propagates to the right and undergoes evanescent 
decay to the left, (c) The time-averaged electric field intensity 
along the mid-line of the waveguide (dashed line in b). 



In this Letter, we show that one-way modes can be gen- 
eralized to photonic crystals with gyrotropic constituents 
without the restriction of having Dirac points in the 
bandstructure or the use of gyroelectric materials. To 
do this, we derive an analytical mapping between the 
electromagnetic modes in the MO photonic crystal and 
the wavefunctions of a nonrelativistic electron in a QH 
system. This links photonic one-way modes to QH edge 
states via the Hatsugai condition, which relates the num- 
ber of edge modes in a band gap to the sum of Chern 
numbers for all bands below it |9j. Using gyromagnetic 
materials, we design an experimentally-feasible one-way 
waveguide using a 2d MO square-lattice Yttrium-Iron- 
Garnet (YIG) photonic crystal (which lacks Dirac points) 
operating at microwave frequencies. Such a one-way 
waveguide consists of an interface between an MO pho- 
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tonic crystal and a gapped material, such as a regular 
2d photonic crystal, as shown in Fig. Q] With realis- 
tic material parameters, the one-way modes are later- 
ally confined to a few lattice constants, and occupy a 
broad (~ 10%) band-gap with negligible material loss. 
The MO photonic crystal exhibits strong (order unity) 
T-breaking due to the use of gyromagnetic ferrimagnets. 
Using both frequency and time-domain numerical simu- 
lations, we demonstrate that these edge modes are im- 
mune to backscattering. The 2d waveguide structure can 
be mapped into an equivalent 3d structure that supports 
identical field distributions and transport properties [l(| • 
For definiteness, consider an MO crystal consisting of 
a square lattice of YIG rods (e = 15eo) of radius 0.11a 
in air, where a is the lattice constant. An external DC 
magnetic field applied in the out-of-plane (z) direction 
induces strong gyromagnetic anisotropy, with the per- 
meability tensor taking the form 



\i in 
-in /i 
/i„ 



(1) 



With a 1600 Gauss applied field, the tensor elements in 
YIG at 4.28 GHz are k = 12.4/i and /i = 14/x [H|- 
For now, we neglect the effect of material dispersion and 
loss, assuming a frequency-independent permeability ten- 
sor with real- valued k and \i. We will later show that our 
results are not substantially affected by these effects. 

Due to the presence of magnetic anisotropy, we must 
adapt the conventional band theory of photonic crystals 
Q to this system: we eliminate the magnetic field from 
Maxwell's equations to obtain the master equation 



V x (/u- x (r)V x E) = e(r)w 2 E. 



(2) 



Here, the inverse permeability tensor and the scalar 
permittivity e are both functions of position, and u> is the 
mode frequency. This equation can be cast in Hamilto- 
nian form if we define the Hermitian inner product as 



(Ei|E 2 ) = / d 2 r e(r)EJ • E 2 . 



(3) 



The existence of one-way edge modes relies critically 
on T-breaking in MO crystals and its effects on the topo- 
logical properties of the Bloch bands. This effect is char- 
acterized by the Chern number, a quantity that has been 
extensively studied in QH systems [H, Ej]. The Chern 
number of the n-th photonic band is 
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where the fc-space integral in ([4]) is performed over the 
first Brillouin zone, and the Bloch-function inner prod- 
uct in ([5]) is defined similarly to ([3]) with the integral per- 
formed over the unit cell. The key properties of the Chern 
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FIG. 2: (Color) Construction of a MO photonic crystal sup- 
porting one-way edge modes. The crystal consists of a square 
lattice of YIG rods (inset in b, with e — 15eo and r = 0.11a) in 
air. (a) Band diagram with zero DC magnetic field (/i = /Uo, 
k — 0). The relevant quadratic degeneracy point is indicated, 
(b) Band diagram with a 1600 Gauss +z DC magnetic field 
(fj, — 14/io, K = 12.4/io)- The degeneracies are lifted, resulting 
in the given non-zero Chern numbers (red numbers) . (c) Con- 
tour plot for the second band. Although the MO crystal has 
only Ca symmetry, the band structure retains an accidental 
Ci v symmetry, with an irreducible Brillouin zone identical to 
that of a non-MO crystal. The Chern numbers are calculated 
by integrating Eq. [5] along the boundary of the first Brillouin 
zone, (d) The corrected band gap with material dispersion. 



number are that (i) it is always an integer, (ii) the sum 
of the Chern numbers over all bands is zero, and (iii) the 
Chern number of every band is zero if the Hamiltonian is 
T-symmetric Property (i) implies that as one adi- 
abatically tunes the Hamiltonian (e.g., varying the per- 
meability tensor by tuning the DC magnetic field bias), 
the Chern number of a given band changes, if and when 
it does so, abruptly. In fact, this abrupt change occurs 
at the critical point where the band becomes degenerate 
with a neighboring band, at discrete "degeneracy points" 
in fc-space. When the system is tuned past this critical 
point, the degeneracies are lifted; a bandgap opens and 
the band's Chern number changes by an integer AC (In 
other words, the relevant bands acquire nonzero Chern 
numbers as soon as a nonzero bias is applied. However, 
a large bias is necessary to support a wide bandgap and 
a tightly-confined edge mode.) Due to property (ii), a 
Chern-number change in one band is accompanied by an 
equal and opposite change, —AC, in the Chern number 
of the neighboring band. Hatsugai has shown that in 
a lattice QH system, this process creates a total of AC 
different edge states at the system's boundary [§]. 

This result, together with property (iii), implies a sim- 
ple and general strategy for constructing photonic crys- 
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tals supporting one-way modes. We begin by searching 
for a T-symmetric band structure with a pair of pho- 
tonic bands degenerate at discrete fc-points; the require- 
ment for Dirac points where the bands touch in a linear 
fashion, which was discussed by Raghu and Haldane Q, 
is not strictly necessary. For example, Fig. [5] shows the 
TM band diagram for a square lattice of YIG rods in air, 
containing a quadratic M-point degeneracy between the 
second and third bands. This degeneracy is lifted upon 
applying an external DC magnetic field, which introduces 
magnetic anisotropy in the high-index rods and breaks T 
(as well as inversion and mirror symmetries.) The sep- 
arated bands acquire non-zero Chern numbers, causing 
one-way modes to appear at the edges of the sample. 

The Chern number argument for the existence of pho- 
tonic edge modes depends on a crucial assumption: that 
Hatsugai's relation between edge states and Chern num- 
bers [9j, which was derived using a lattice QH model, 
applies to the photonic system. Although this has not 
been formally proven, there is an interesting relationship 
that allows us to map the photon states in gyromagnetic 
crystals to electron wavefunctions in a family of QH sys- 
tems. To see this, let us return to ([2]) and write /i. _1 (r) 
in terms of its component scalar functions: 
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In terms of the component functions in ([Tj). fj, 



u- 1 = 

The equations for E x 



\ij (/j, 2 — k 2 ) and rj = — «/ (fi 2 — n 2 ) 
and E y decouple from E z ; for TM states (E x = E y = 0) , 

[- V 2 + (V In A - ijiz X Vt?) • V - fleuj 2 ] E z = 0. (7) 

Expressing this in terms of ip = E z y/p,, we obtain 



V + iA(r) +V(r) 



i> = o, 



where 



A = - z x Vr/, 



(8) 
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V = i (|VlnA| 2 + l/iVr?! 2 ) -iv 2 lnA-M^ 2 -(10) 

This is the equation for zero-energy wavefunctions of a 
non-relativistic particle in periodic vector and scalar po- 
tentials A(r) and V(r). Increasing to corresponds to in- 
creasing the depth of the scalar potential well in the third 
term of (|10|) . shifting the spectrum downwards relative 
to the zero of energy. With suitable boundary condi- 
tions, this mapping holds for both edge modes and bulk 
modes. Thus, for each value of uj, the existence of un- 
paired edge modes, as well as their spatial characteristics, 
can be mapped to a similar problem in a QH system. 

To complete the single-mode one-way waveguide, we 
interface the MO crystal with a upper cladding that sup- 
ports no bulk modes at the frequency of the second MO 



crystal bandgap. This can be achieved using a metal 
slab, or using a regular photonic crystal with an aligned 
bandgap. The topological nature of the edge modes @ 
ensures that the interface structure is unimportant, as 
long as it is sufficiently narrow that higher order modes 
are avoided. Here, we use a square lattice of high-index 
alumina rods (r = 0.106a, e = 10eo) in air, tilted 45 de- 
grees to match the bandgap frequency (Fig. Q})). Note 
that the regular crystal has zero Chern numbers for all its 
photonic bands, so there is only one forward-propagating 
mode at the interface. With a 1600 Gauss bias, we obtain 
a wide MO bandgap of 10% relative size and very low dis- 
persion, with the edge mode confined to a few lattice con- 
stants. Using a finite-element frequency-domain scheme 
and a steady-state excitation, we observe an evanescent 
mode decaying exponentially along the backward direc- 
tion of the waveguide (Fig. [TJ:). 

One of the unusual properties of the one-way edge 
modes is the complete suppression of backscattering, 
which is an optical analog of the dissipation-free trans- 
port of edge electrons in QH systems. Numerical simu- 
lations show that these edge modes are immune to scat- 
tering from extremely large defects. For instance, Fig. [3] 
shows the results of simulations with a slab of perfect 
electrical conductor (PEC) of width 3a and thickness 0.2a 
inserted into the waveguide. In a conventional waveguide, 
such a drastic defect would almost completely block the 
guided mode. In the one-way waveguide, a steady-state 
source operating at the mid-gap frequency 0.555 • (2irc/a) 
excites a one-way mode that circumvents the PEC de- 
fect, with 100% power transmission throughout the MO 
bandgap. This happens because the defect creates a new 
interface waveguide between the PEC and the MO crys- 
tal. Thus, it only alters the phase response, which is 
partly due to the delay incurred by traversing the length- 
ened interface. We corroborated this result with a time- 
domain calculation in which a temporal Gaussian pulse 
with a spectral bandwidth of 50% of the bandgap and car- 
rier frequency 0.555 • (27rc/a) is launched into the waveg- 
uide. Regardless of the presence of the defect, the pulse 
passes through the waveguide with no perceivable change 
in amplitude or pulse width; since the one-way mode has 
approximately linear dispersion relation at mid-gap fre- 
quencies (Fig. [TJ, sharp corners do not contribute sig- 
nificantly to chromatic dispersion. The increased transit 
time is in agreement with the change in group delay. Al- 
though embedded sources were used in these simulations, 
we have verified in a separate set of simulations that it 
is also easy to couple to the one-way mode using waves 
incident on the system boundary. 

While our analysis has thus far been limited to 2d 
structures at microwave frequencies, similar one-way 
waveguides can be realized in practical 3d structures and 
at higher frequencies. In the 2d analysis, the dielec- 
tric structure and electromagnetic field extend uniformly 
along the z direction. In a 3d structure, TM modes can 
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FIG. 3: (Color) Back-scattering suppression in the one-way 
waveguide. Inserting a slab of perfect electrical conductor 
with thickness 0.2a (black rectangle) causes the propagating 
modes to circumvent the defect, maintaining complete trans- 
mission. E z is plotted for (a) t — 0, (b) t — 0.25Tn, and (c) 
t — 0.5Tb, where To is the optical period, (d) Time-domain 
simulation results for a temporal Gaussian pulse with spec- 
trum contained in the bandgap. The electric field amplitude 
is plotted at the source point (black), at the same transverse 
position 13 lattice constants downstream along the waveguide 
in the absence of a defect (red) , and at the same point with an 
intervening defect (blue). The pulse is completely transmit- 
ted regardless of the existence of the defect, (e) Transmission 
and phase shift plots for the time-domain simulation. 



be truncated by introducing a PEC (for which a metal 
slab is a good approximation at microwave frequencies) 
in the x-y plane, without affecting the field distribution. 
An array of finite high YIG rods and alumina rods be- 
tween two metal slabs supports TEM modes with field 
distributions identical to the TM modes in our 2d sys- 
tem 

[13, EH- Finally although gyromagnetic effects are 
generally limited to the GHz range by the availability 
of high-field magnetic bias, artificial magnetic resonance, 
such as those exploited in metamaterials, might be incor- 
porated to synthesize the gyrotropic responses at THz or 



ity in a real microwave ferrite. With a gyromagnetic 
linewidth of 0.3 Oe and a dieletric loss tangent of 0.0002, 
typical in commerically available monocrystalline YIG 
[17j |. the complex propagation constant is (—0.359 + 
0.0001i)(27r/a) in the case of Fig. 3. The imaginary part 
corresponds to a decay length of 300 lattice constants, 
far exceeding practical structural dimensions. We have 
also neglected the frequency dependency of the gyromag- 
netic permeability. When correcting this approximation 
with a complete gyrotropic lineshape [ll| , we find the MO 
bandgap reduces to around 6%. Although such reduction 
narrows the operational bandwidth, we find the disper- 
sion has no impact on the back-scattering suppresion or 
the confinement of the edge mode. We envision many ap- 
plications for these one-way waveguides, including defect- 
tolerant slow-light systems Q, microwave isolators, high- 
Q channel add/drop filters, and all-pass filters. It might 
also be possible to further exploit the analogy between 
this system and the QH effect, by for instance reproduc- 
ing the fractal "Hofstadter butterfly" spectrum [n| in 
the bulk modes of the MO photonic crystal. 
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even infrared frequencies [lj, 13, 13] • 

The above results are not substantially altered by 
the material losses and frequency-dependent permeabil- 
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